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p) • Lepton number violation (LNV) can be induced by Majorana neutrinos in four-body decays of 

I) : 

. the neutral B meson and the top quark. We study the effects of Majorana neutrinos in these 

, l^-^l = 2 decays in an scenario where a single heavy neutrino can enhance the amplitude via 
> ■ 

Q>^ ■ the resonant mechanism. Using current bounds on heavy neutrino mixings, the most optimistic 
O i 

. branching ratios turn out to be at the level of 10 ® for — ?• D^e e vr^ and t — )• bl~^l~^W 



od 
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X 



decays. Searches for these LNV decays at future facilities can provide complementary constraints 
on masses and mixings of Majorana neutrinos. 
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I. INTRODUCTION 



After being established that neutrinos are massive and mixed particles [l-^, one of the 
most interesting current issues in flavor physics is to elucidate if neutrinos are Dirac or Ma- 
jorana fermions j4|. If neutrinos turn out to be Majorana particles, important consequences 
as lepton number- violating (LNV) processes js], Q and further sources of CP violation be- 
come possible 

Q-|lo|- Searches for LNV processes (where the lepton number is violated in 
two units, AL = 2) in dedicated low energy experiments as neutrinoless double beta decays, 
have led to very strong constraints on the effective mass of light Majorana neutrinos llMlS] 
since the rates for these processes are driven by the effective mass parameter {rriee) isl- On 
the other hand, very restrictive bounds on (rnw) can be obtained by combining neutrino 



oscillation data [17] , cosmological bounds 18l-l2l| and tritium beta decay 22|. Interestingly, 
these sub-eV bounds on the scale of effective Majorana masses are at the sensitivity reaches 
of current experimental projects 



23 



26|. 



As it has been extensively discussed by many authors, the existence of very light neutrinos 
may find a natural explanation by means of heavy neutrinos via the see-saw mechanism 
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33|. Heavy neutrinos naturally appears in some extensions of the Standard Model and 
may play an important role in cosmology and various particle physics and astrophysical 



processes 3J, ISSj. The possibility to observe the effects of heavier neutrinos, accessible 



in the kinematical range of current experiments, is very exciting as they can induce large 
rates for AL = 2 decays through the mechanism of resonant enhancement 35|]. Indeed, 
the appearance of sterile neutrinos with masses in the range of hundreds of MeV's to a 
few GeV's is possible in scenarios of dynamical electroweak symmetry breaking as shown for 



instance in 



36N38| . By means of the resonant mechanism, neutrinos with these intermediate 



mass scales can produce an enhancement in the three-body AL = 2 decays of pseudoscalar 
mesons — )■ l^l^M2 and the tau lepton r~ — )■ l'^ M2 '■, these decay processes have 



been extensively studied by many authors 35|, l39N43| in the cases where final state hadrons 



can be pseudoscalar or vector mesons. So far, some experimental uppers bounds have been 
reported in refs. 



44 



47j; very recently, by using 36 pb— 1 of integrated luminosity, the 



LHCb collaboration has been able to derive upper limits for LNV charged B meson decays 
5(5+ i^~(7r~)/i+/i+) < 5.8(5.4) X 10"^ |48[. These studies are expected to be extended 
by the LHCB experiment by including the — )• D^^^ji'^n'^ ^ Z)°/i+/x+7r^ decay modes |48| . 



which toghether with similar analyses that can be performed at the SuperB Flavor Factories 



491 makes very attractive the studies of LNV B meson decays. Similarly, like-sign dileptons 



may be produced via the resonance enhancement mechanism in four-body decays of top 
quarks and W gauge bosons, as it has been investigated for instance in Refs. 50l-l52|. 

In the present paper we consider the four-body decays of neutral B mesons, B^ — )■ 
D^l^l^n^ with I = e, /i, in the favored scenario of resonant neutrino enhancement. The 
dynamics of this four-body decay involves the transition B ^ D form factors and is different 
from the one driving the three-body decays of mesons and tau leptons which involve the 
meson decay constants. To our knowledge, these AL = 2 decays of neutral B mesons have 
not been investigated before neither from a theoretical nor from an experimental point of 
view. In addition, we also consider and update the analogous four-body t — ?■ bl~^l~^W~ 
decays (/ = e, fi, r), which was previously studied in [50], since one naively expects it 
can be largely enhanced due to the resonances in the virtual W boson and heavy neutrino 
exchanges. 



II. FOUR-BODY AL = 2 DECAYS OF HEAVY FLAVORS 

The Feynman diagrams that describe the LNV decays of the top quark and neutral B 
meson are shown in Figures 1 and 2, respectively. Contributions containing the properly 
antisymmetrized contributions due to the exchange of identical leptons in the final state 
must be added to these diagrams. 



Following previous studies j35|,|50|, we consider a model with three left-handed SU(2) lep- 
ton doublets = (i/q, la)L, (a = 1, 2, 3), and n right-handed singlets NbR (6 = 1, 2, ■ ■ -n). 
n the basis of mass eigenstates, the charged current interactions of leptons are given by 



35|: 



Cf^ = -jjf^t I E E ^^rnVrnl^Pd + E E UlmN^^Pd ) + h. C. (l) 



l=e 771=1 l=e m=l 



where Pl = (1 —75)72 is the left-handed chirality operator, g is the SU{2)l gauge coupling, 
ifj'^ = Cip'^ is the charge conjugated spinor, and Vim iUim) denotes the light (heavy) neu- 
trino mixings; the subscript m refers to the mass eigenstate basis entering the diagonalized 
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FIG. 1. Feynman graph for the t — )• bl^l^W decay. 
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FIG. 2. Feynman graph for the AL = 2 neutral B meson decay. 



Majorana mass term for neutrinos j35| : 
C = - 



\m=l m=4 / 



(2) 



In the phenomenological appUcations of the present paper, we will assume that only one 
heavy neutrino with mass rriN and charged current couplings Um to leptons, dominates the 
decay amplitudes via the resonant enhancement mechanism. 

The kinematics of four-body decays can be described in terms of five independent vari- 
ables. In our convention of momenta and masses they are defined as 

with p^ = and pj = mf. We choose the set of independent variables as {si2, S34, 9i, 9^, 0} 
which have the following geometrical meaning 53|] (see Figure 3): 



z 




y 



FIG. 3. Kinematics of four-body decays in the rest frame of the decaying particle, Yli=iPi = 0. 
We have defined pij = Pi+ Pj, such that pi2 + P34 = 0. 

• S12 = (pi +^2)^, is the invariant-mass of particles 1 and 2; 

• S34 = (p3 +P4)^, is the invariant-mass of particles 3 and 4; 

• 61 {O3), is the angle between the three- momentum of particle 1 (particle 3) with respect 
to the direction of pi2 = pi+p2 (respectively, P34 = + p^) defined in the rest frame 
of the decaying particle; 

• is the angle between the planes defined by particles (1,2) and (3, 4) also in the rest 
frame of the decaying particle. 

With this choice of kinematics, the differential decay rate in the rest frame of the decaying 
particle can be written as: 

= TTprin^l-^l'^ ■ -^dsi2ds34dcos9idcos9sd(f) , (3) 
4(47rj°M'^ n\ 

where P12 {P34) is the velocity of particle 1 (particle 3) in the center of mass frame of particles 
1 and 2 (3 and 4) and X = X^/^(M^, su, S34) /2, with X{x, y, z) = x^+y'^+z'^-2xy-2xz-2yz. 
Finally, |A1p is the spin-averaged squared amplitude of the four-body decay. In our case of 
two identical leptons in the final state, n = 2. 
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III. B DUtt decays 



Let us first consider the B^{p) — > {p2)l~ {p^)t^^ iP'i) decays (see Figure 2), where 

Pi denote the four- moment a of final state particles. In the range of neutrino masses ttin 
where the resonance effects dominate the decay amplitude, the diagrams of Figure 2(c) 
and 2(d) will give very small contributions. In addition, we note that diagram 2(b) is 
suppressed with respect to 2(a) due to smaller Cabibbo-Kobayashi-Maskawa (CKM) factors 
{\VubVcd/ {VcbVudjl ~ 0.02). Therefore, we keep the diagram shown in 2(a) as the dominant 
contribution. 

The properly antisymmetrized decay amplitude is given by: 

M = GlV,,VUD{vi)\crh\B{p)) ■ u{v2) [Vn{V2)i,1u + Vn{Vz)i,i,] u%ps) {lUvX) (4) 

where Vij denotes the ij entry of the CKM quark mixing matrix, Gj? is the Fermi constant 
and = 130.4 MeV is the 7r+ decay constant. 

In the above expression we have defined the factor 

-PNipd = ^ , (5) 

[Q - Pi) + trriNTN 

where Q = p — Pi = P2 + Ps + P4 is the momentum transfer and Uuy denotes the heavy 
neutrino mixing defined in Section 2. F^v represents the decay width of the heavy neutrino 
which depends on the decay channels that can be opened at the mass mAr; it allows to keep 
finite the amplitude when {Q — Pi)"^ = mj^. As it was pointed out above, in this paper we 
will assume that only one heavy neutrino falls in the resonance region of the B meson and 
top quark decays, thus it will give the dominant contribution to the decay amplitude. The 
mixings of the heavy neutrino with the three charged leptons will be taken as the currently 



most restrictive bounds as reported in Ref. 



54| 



Setl: |?7eivP < 3 X 10"^ |f/^jvP < 3 x 10~^ |f/^^p < 6 x lO'^ . (6) 

The hadronic matrix element in Eq. (4) is given by: 

{D^p,)\-n,bmp)) = ((p + pi), - ^^p^ g,) F,{t) + ^^;^^ Q,Fo(t) , (7) 

where t = Q^. For the purposes of a numerical evaluation, we will use two common 
parametrizations of the form factors Fi^o{t), namely the one provided by the Wirbel-Stech- 
Bauer (WSB) model 5^: 

T7WSB/.\ _ 
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FIG. 4. Neutrino decay width for neutrino masses relevant to produce resonant enhancement in 
B° D+1-1-TT+ decays. 



where F^sb^q) = F^^^{0) = 0.69, mi- = 6.34 GeV and mo+ = 6.8 GeV |55| and, just 
for comparison, we will use also the parametrization provided by the covariant light front 
(CLF) model |56|: 

F^'^'iO) pCLF..^ F,^^^iO) 



1 -ai(t/m|) +6i(t/m|) 



2 •\2 ' 



ao{t/m%) + bo{t/m%) 



2 \2 ' 



(9) 



where FfLF(o) = ^^^^^(O) = 0.67, = 1.25, h = 0.39, = 0.65 and bo = 0.0 |56j. 

The decay width of the intermediate neutrino state is obtained by adding up the contri- 
butions of all the neutrino decay channels that can be opened at the mass ttin |35| : 



(10) 



where ruf. in the argument of the step function are the masses of the final state particles 
in the neutrino decay channel /. The dominant decay modes of the neutrino in the range 
of masses relevant for resonant B meson decays are the following: PP^, uiP^, PV^, viV^, 
Iflfui^, z/^j/^/J, and vi^i'i'., where Z, li, h = e, /i, r, and P {V) denotes a pseudoscalar 
(vector) meson state. The expressions for the partial decay rates of these channels can be 
found in Appendix C of Ref. 35[. We have re-evaluated the decay width of the neutrino 
which is plotted in Figure 4 for neutrino masses m^v in the range where it can produce a 
resonant enhancement of the B meson decay amplitude. The decay width Fjv is so tiny that 
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FIG. 5. Branching ratio of — t- D^l I vr^ decay as a function of mj\[. The dashed (sohd) hne 
corresponds to the electronic (muonic) channel. 



the Narrow Width Approximation (NWA) of Eq. (5), 



lim VN{Pi) = -inniNUi^S ({Q - Pif - m%) 

Tat— *-0 



(11) 



is required to perform the numerical integrations to compute the decay rates from Eq. (3). 

In Figure 5 we plot the branching ratios of — t- D^l^l^n^ decays for the electron 
(dashed-line) and muon (solid-line) channels as a function of the neutrino mass rriN- These 
plots were obtained by using the WSB model 55|] for the form factors Fi o(t). The branching 
fractions reach their maximum values for neutrino masses that are close to the threshold for 
the mi + m,r production and they decrease for increasing values of mAr. 

In Table I we show the largest possible values of the branching ratios of B'^ — ^ D~^1~1~tt~^ 
decays (Z = e, /i), which correspond to the lower range of neutrino masses. We have 
evaluated these results for the two different form factor models mentioned in Eqs. (8) and 
(9). Although the predictions for the form factors in the two models exhibit large differences 
in the full range of the momentum transfer t, the integrated rates differ only at the level 
of 30% for almost all values of neutrino masses, both in the electronic and muonic decay 
channels. The largest possible values of the branching fractions shown in Table I are of 
the same order as the ones corresponding branching ratios reported for the B^ — t- P^l^l 
decays 



4J] and their study can be useful to get further constraints on the heavy neutrino 

8 



' '^N V V ) 


VV 0-L> 








' ' \ V ) 


VV 0-L> 








i ( u 


Z.D X lU 


-6 


0.4 X lU 


-6 


zou 


O.U X iU 


-7 


o.y X iU 


-7 


190 


2.8x10" 


-6 


3.6x10" 


-6 


270 


4.1x10" 


^7 


5.4x10" 


-7 


200 


2.6x10" 


-6 


3.4x10" 


-6 


300 


3.4x10" 


-7 


4.3x10" 


-7 


220 


1.5x10" 


-6 


2.0x10" 


-6 


400 


1.4x10" 


-7 


1.9x10" 


-7 


250 


7.3x10" 


-7 


9.7x10" 


-7 


500 


7.0x10" 


-8 


1.0x10" 


-7 


300 


2.5x10" 


-7 


3.3x10" 


-7 


600 


4.0x10" 


-8 


6.0x10" 


-8 



TABLE I. Branching ratios for — t- D^l £ vr"'" decays using the Set I of the heavy neutrino 
mixings. WSB 



55l | and CLF [5S] refer to the form factor models for the B ^ D transition. 



mixings. 



IV. LNV TOP QUARK DECAY 

LNV transitions with AL = 2 has been studied also at higher energies. The t — > hlfl~^W~ 
decay (and its crossed — )■ If^'j^ 2 jets channel) has been considered previously in Ref. 



50|; similar top quark decays that also include the final W gauge boson decay into two 



jets were studied in 511 • The top decay can be resonantly enhanced if the heavy neutrino 
mass lies in the range mw + mi < < frit — rrii, — mi. In addition (see Figure 1) we 
can expect an enhancement of the top quark decay amplitude due to the virtual W boson 
emitted from the to p qu ark vertex which can be produced also in a resonant way. As it was 



emphasized in ref. [50|, the final state W boson in this AL = 2 top quark decay has the 
'wrong' charge signature when compared to the dominant t — )■ bW^ decay. In this Section 
we provide an update for this same-sign dilepton decay channel in top quark decay by using 
the neutrino mixings provided in Eq. (6); furthermore, our results for this decay channel 
provide a test for the particular kinematics that we use in our calculations and that was 
described in Section 2. 



For the purposes of comparison with previous results on AL = 2 top quark decays 50|, 



we will evaluate our results using, in addition to the Set I of values given in Eq. (6), the 



following set of neutrino couplings 57 1 



Setll: \UeN\ < 12 xlQ--", |[/^^|^ < 9.6 x 10"^ If/^^vf < 16 x lO^'^ . (12) 

The decay amplitude for t{p) — )■ b{pi)l'^ {p2)l^ {p3)W~ (pi) corresponding to the diagram 
of Figure 1 is given by: 



V2 \V2, 

xm(p2) [VNiP2Hl. + P7V(P3)7.7/.] u'iPs) " (p^) , (13) 

where D^{Q) = ii-gpf" + QPQf'/m'^)/{Q'^ - + imw^w), with Q = p ~ p^, denotes the 
resonant W boson propagator in the unitary gauge, e'^{p4) is the polarization four- vector of 
the W~ boson and VNiPi) was defined in Eq. (5). 

The total decay width of the neutrino for the range of neutrino masses giving rise to 
the resonance enhancement, is determined from the following set of two-body final states: 

— l^W^ , uiZ^ and uiH. The expressions for the total width by neglecting the charged 



lepton masses is given by [35|, |50[ : 



= [2im% + 2ml,)Xw + {m% + 2m|)Xz + m^Xn] , (14) 



where Xi = (m^ — mf)'^6{mN — rrii) for i = W^, Z and H bosons. As long as the 
neutrino mass increases, the total decay width r^r also grows because of the neutrino mass 
dependence and also because new decay channels are opened. For neutrino masses relevant 
for top quark decays, Fjv is several orders of magnitude larger than for B meson decays 
and a straightforward evaluation of the five-dimensional integration of Eq. (3) can be done 
without numerical complications; at the same time, the neutrino width is small enough that 
it allows also the use of the NWA approximation, Eq. (11), to integrate the phase space. 

In Figure 6 we plot the branching fraction (normalized to It/iArl"^) of t — )■ bl^l^W~ decays 
as a function of m^v. The phase space and the squared amplitudes are almost insensitive to 
the masses of different leptonic channels in the final state. We use as inputs: rrit = 172.0 



GeV, tuh = 120 GeV, the leptons and gauge bosons masses given in ref. [4J], and the 
neutrino mixings given in Eq. (6). The solid-line represents the branching ratio that is 
obtained from the five- dimensional integration of Eq. (3). The dashed-line is the result 
obtained by using the NWA method. Clearly, the results obtained by using these two 
methods are almost identical except for very small differences appearing at the upper values 
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FIG. 6. Normalized branching ratio of t — t- bl^l^W~ decay as a function of m^v- The sohd-Hne is 
obtained by a straightforward integration of Eq. (3), and the dashed-Une corresponds to the NWA 
for the neutrino width. 

of neutrino masses that can produce the resonant enhancement. As it was pointed out above, 
the results shown in Figure 6 were obtained by using the Set I of neutrino mixings given 
in Eq. (6). Resuhs for a new set of neutrino mixings can be obtained by muhiplying the 



results shown in Figure 6 by T/(Xli 1^/ 



/"new 1 2 
IN I ) 



In Table II we show the branching ratios of t — > bl^l^W^ {I = e, /i, r) decays for a 
few values of neutrino masses such that the rates have their largest values. For this specific 
range of neutrino masses, the branching ratios turn out to be of order 10^^ ~ 10^^. For a 
given Set of neutrino mixings, the results for different leptonic channels differ basically by 
the rescaling of their fourth power of Uin, as it should be. Just for comparison, we have 
computed the branching ratios by using the Set II of neutrino mixing parameters; our results 
are compared with those in Ref. jsoj] (shown within parenthesis in Table II). Our results 



and those of reference 



50| have similar values for the electronic and muonic channels, but 



they differ in the rr channel by about 30%. Finally, let us comment that we have evaluated 
the branching ratios of t — )■ bl~^l~^W~ for a wider range of the heavy neutrino mass. The 
normalized branching ratio B{t — ?■ bllW)/\Uii\i\^ plotted in Figure 6 drops from 1.7 x 10~^ 
to 8.3 X 10^^^ when the neutrino mass spans from 200 GeV to 2 TeV. 

We close this section with an estimate of the expected sensitivities to the signals of our 
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Set I 


TUN (GeV) 


ee 


MM 


TT 


90 


0.29 


0.29 


1.12 


100 


0.12 


0.12 


0.47 


110 


0.05 


0.05 


0.19 


Set II 


mjv (GeV) 


ee 


MM 


TT 



90 1.48 (1.4) 0.95 (1.1) 2.55 (1.9) 
100 0.62 (0.6) 0.40 (0.5) 1.08 (0.8) 



5C| 



TABLE II. Branching ratios (in 10 ^ units) for t bi^£'^W~ decays. Results of Ref. 
corresponding to the Set II of neutrino mixings are shown within parenthesis. 



AL = 2 decays. Using as a (conservative) reference the ~ 450 x 10^ BB pairs accumulated 
by the BABAR collaboration at the T(4S'), we can provide an estimate of the sensitivity to 
the branching ratio of B^ — > D^l^l^n^ decay at the B factories by using the K'tt^tc'^ mode 
to reconstruct the charged D meson. By assuming a 70% efficiency for the identification 
and reconstruction of each of the six charged tracks in the final state one can reach a 
sensitivity of ~ 2.0 x 10^^ which can test some range of our upper limits for the di-leptonic 
channels. Of course, this optimistic estimate does not include the combinatorial background 
for this detection channel, although it can motivate a more detailed study of backgrounds 
and efficiencies (note also that, under similar assumptions, slightly better sensitivities can be 
reached using Belle data which is about the double of B meson pairs produced by BABAR). 
Improved sensitivities can be obtained at the Super Flavor Factories |49| which are expected 
to accumulate a larger data set by a factor of 50 to 75 with respect to B factories. Also, as 
it was mentioned in the Introduction, there are good perspectives at the LHCb Experiment 
which can provide sensitive improvements on these LNV decays in the dimuon channel 



based on the analysis already done in the case of LNV searches in B^ meson decays 48 1. 
Regarding the top quark AL = 2 decays, the sensitivies that can be reached at the Large 
Hadron Collider (LHC) are not sufficient to test our predictions, as it has been discussed in 

12 



Ref. |50|- Only in the eventual case of an upgraded Super-LHC Experiment, which should 
increase the LHC luminosity by a factor of 10 [58], one can expect that branching ratios of 
10^^ for t — )■ bl^l'^W' decays could be accessible. 



V. CONCLUSIONS 

The existence of lepton number- violating transitions with AL = 2, is considered to be the 
cleanest manifestation of neutrinos as Majorana particles. Direct searches of these decays 
in neutrinoless double beta nuclear decays indicate that masses for very light Majorana 
neutrinos lie at the sub-eV scale. Similarly, in the framework where only three light neutrinos 
exist, the Majorana masses of all neutrinos are strongly constrained from current oscillation, 
cosmological bounds and tritium beta decay. 

In the present paper we have studied the potential of heavy flavor four-body AL = 2 
decays to shed some light on the masses and couplings of heavier Majorana neutrinos. If 
the masses of such neutrinos produce a resonance enhancement of these heavy flavor decays, 
the corresponding branching ratios turns out to be, in the most optimistic cases, at the level 
of 10~^ for neutral B meson and top quarks decays if the most restrictive current bounds 
on neutrino mixings (Set I) are used. These branching fractions of four-body neutral B 
meson decays are at the level of the upper bounds obtained in experimental studies of B^ 
three-body decays. Thus, their searches at current and future experimental facilities can be 
helpful to provide complementary constrains to the ones derived from three-body decays of 
r leptons and charged B and D mesons. 
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